Universal Quantum Computation with the z/ = 5/2 Fractional Quantum Hall State 
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We consider topological quantum computation (TQC) with a particular class of anyons that are 
believed to exist in the Fractional Quantum Hall Effect state at Landau level filling fraction v — 5/2. 
Since the braid group representation describing statistics of these anyons is not computationally 
universal, one cannot directly apply the standard TQC technique. We propose to use very noisy 
non-topological operations such as direct short-range interaction between anyons to simulate a 
universal set of gates. Assuming that all TQC operations are implemented perfectly, we prove 
that the threshold error rate for non-topological operations is above 14%. The total number of 
non-topological computational elements that one needs to simulate a quantum circuit with L gates 
scales as L(logL)'^. 



I. INTRODUCTION 

One of the most important results in the theory of 
fault-tolerant quantum computation is the threshold the- 
orem. It asserts that ideal quantum circuits can be effi- 
ciently simulated by noisy circuits if an error rate of indi- 
vidual gates is smaller than a certain constant threshold 
value S, see [1-4]. Estimates of S vary from 10"'^ - 10^'' 
for a local architecture [5, 6] to 10~^ — 10~^ for non-local 
gates [4, 7]. With the present technology these rates are 
hardly achievable by any real device. Moreover, for prac- 
tical computations it is desirable to have an error rate 
much smaller than 5, otherwise one may need too many 
concatenation levels and the simulation overhead may 
become too large. 

These challenges can be overcomed (at least partially) 
in the Topological Quantum Computation (TQC) scheme 
developed by Kitaev, Freedman, and Preskill [8-10]. It 
makes use of the fact that elementary excitations of 
some two-dimensional many-body quantum systems are 
anyons — spacially localized quasiparticles with unusual 
exchange statistics described by non-trivial representa- 
tions of the braid group. For the purposes of TQC 
one needs non-Ahelian anyons (corresponding to multi- 
dimensional braid group representations). A computa- 
tion is carried out by creating pairs of anyons from the 
ground state, separating them far apart, transporting in- 
dividual anyons adiabatically around each other, and fi- 
nally fusing pairs of anyons together. A list of particle 
types produced in the fusion is the classical outcome of 
the computation. An error rate of individual gates in 
TQC is expected to be much smaller than 5. 

A physical system that may serve as a platform for 
TQC is a two-dimensional electron gas in the Frac- 
tional Quantum Hall Effect (FQHE) regime. The FQHE 
plateau at the filling fraction v = hjl was observed by 
Willett et al. [11] in late eighties. Shortly after that 
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Moore and Read [12] developed a theory predicting that 
elementary excitations of the v = 5/2 state are non- 
Abelian anyons. The corresponding braid group repre- 
sentation was found by Nayak and Wilczek [13]. For the 
sake of brevity we shall refer to the anyons existing in the 
V = b/2 state as Ising anyons (their exchange statistics 
can be described by monodromy of holomorphic correla- 
tion functions of the 2D Ising model [12]). 

From the experimental point of view, Ising anyons 
have many favorable properties. A large quasiparticles 
gap (estimated as A > 100 mK in [14]) suppresses ther- 
mal creation of 'stray' particles, while non-zero electric 
charge permits control of anyons using electrostatic gates. 
Besides, one can take advantage of the well-developed 
FQHE experimental technology. An experimental setup 
for controlling Ising anyons and testing their statistics 
has been recently proposed by several authors [15-18]. 
An error rate for the one-qubit operation has been 
estimated as lO^'^" in [15]. 

The only fact that prevents one from using Ising anyons 
for TQC is that the braid group representation describ- 
ing their statistics is not computationally universal. We 
shall see that one can easily compute an amplitude of any 
braiding process, see Section HI. Loosely speaking, TQC 
with Ising anyons is an intersection of two computational 
models known to be classically simulatable: quantum cir- 
cuits with Clifford gates [19-21], and Fermionic Linear 
Optics [22-24]. Therefore, Ising anyons offer only reliable 
storage of quantum information and reliable implemen- 
tation of a certain non-universal gate set, see Section III 
for details. 

The goal of the present paper is to argue that this 
drawback is not as serious as it might seem. We show 
that a universal gate set can be simulated by standard 
TQC operations, i.e., adiabatic transport and fusion of 
anyons, and very noisy non-topological operations, such 
as direct short-range interaction of anyons. The latter 
can be thought of as a tunneling process in which two 
anyons exchange a virtual quasiparticle. It can be im- 
plemented by transporting two anyons sufficiently close 
to each other, waiting for an appropriate period of time, 
and then returning the anyons to the original positions. 
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Another example of a non-topological computational el- 
ement is a two-point contact interferometer proposed 
in [16-18]. It has the geometry of a Hall bar with two 
constrictions, such that quasiparticle tunneling occurs 
between two edge currents on the opposite edges of the 
bar. The tunneling current is sensitive to the total topo- 
logical charge of anyons trapped inside the interferometer 
loop. We will show that the short-range interaction and 
the two-point contact interferometer together with TQC 
operations provide a universal gate set. 

Our main result concerns the threshold error rate of 
non-topological operations. To avoid propagation of er- 
rors we apply all non-topological operations before the 
computation itself to prepare a supply of "computation- 
ally universal" ancillary states from the vacuum. In our 
scheme there will be two types of ancillary states: a four- 
particle state 1 04) and an eight-particle state lag)- Prom 
the computational perspective, [04) can be identified with 
a one-qubit state 2~^/^(|0) -|- e"/**|l)) (we represent a 
qubit by four quasiparticles). Analogously, lag) can be 
identified with a two-qubit state 2-^^'^{\0, 0) + \l, 1)). One 
copy of [04) together with TQC operations allows the im- 
plementation of the one-qubit 7r/8-rotation. Analogously, 
one copy of lag) allows the implementation of the CNOT 
gate. Summarizing, universal computation can be car- 
ried out by TQC operations if a supply of states |a4) and 
lag) is available. 

Since non-topological operations are not perfect, in 
practice one can prepare only some very noisy ancillary 
states p4 and ps approximating [04) and \as) up to some 
precision. We characterize this precision by two param- 
eters 

€4 = 1- {ai\p4\a4) and eg = 1 - (ag|pg|ag). 

We prove that the ideal states [04) and |ag) can be dis- 
tilled from many copies of p4 and pg by TQC operations 

provided that 

1. All TQC operations are perfect, 

2. €4 < 0.14, 

3. eg < 0.38. 

A distillation method that we use is a combination of 
"magic states distillation" proposed in [25] and a slightly 
modified version of the entanglement purification proto- 
col of Bennett et al. [26, 27]. 

Summarizing, if one can prepare the states [04), |ag) 
accurately enough, such that the conditions above are 
satisfied, then any quantum computation can be efii- 
ciently simulated by Ising anyons. The overall simulation 
requires only poly-logarithmic overhead. Specifically, the 
number of noisy ancillas p4, ps, and the number of TQC 
operations that one needs to simulate a quantum circuit 
with L gates scales as L (logL)"^. 

In the case when one can meet only the condition 
eg < 0.38, TQC operations allow one to implement any 
Clifford gates (i.e., the CNOT gate, the Hadamard gate, 



and the one-qubit 7r/4- rotation). Though these gates do 
not constitute a universal set, they are sufficient to im- 
plement any error correction scheme based on stabilizer 
codes [28] . Error correction might be needed if one takes 
into account finite error rate of TQC operations (which 
is neglected throughout this paper). 

Our derivation of the threshold conditions 2,3 is based 
on a single assumption regarding the error model char- 
acterizing non-topological operations — they must obey 
the superselection rules of Ising anyons. Accordingly, we 
assume that matrix elements of p4 and ps are non-zero 
only for the vacuum sector (recall that each ancilla is 
prepared from the vacuum). 

The rest of the paper is organized as follows. Section II 
provides the necessary background on Ising anyons. In 
Section III a TQC with Ising anyons is discussed and 
its classical simulatibility is proved. Section IV describes 
a distillation method for the state jag). We show how 
to use ancillas |ag) to implement Clifford group gates in 
Section V. Finally, in Section VI we make use of the 
magic states distillation protocol to simulate universal 
computation. Also the efficiency of the simulation is an- 
alyzed. Some particular non-topological ancilla prepara- 
tion methods are discussed in Section VII. 

II. ISING ANYONS 

A complete specification of any class of anyons is 
rather complicated and involves a lot of data including 
a list of particle types, their fusion and braiding rules, 
S'- matrices, c.t.c, see [29] for a comprehensive review, 
and [16] for a detailed discussion of Ising anyons in the 
context of FQHE. In this section we briefly outline the 
properties of Ising anyons focusing on those relevant for 
quantum computation. 

A. Particle types and fusion rules 

There are two non-trivial particle types in the class 
of Ising anyons. We shall label them by letters a and 
■0. Particles of different type cannot be converted to one 
another (or to the vacuum) by a local operator, thus de- 
scribing superselection sectors of the model. However, 
if one brings two particles close to each other, they can 
fuse into a single one, or annihilate each other forming 
a topologically trivial particle (belonging to the vacuum 
sector). Admissible interconversions of particles are for- 
mally described by fusion rules: 

ip X ip = 1, ijj X a = a, a X a = l-\-il). (1) 

Here 1 stands for the vacuum sector. The most important 

for us is the last rule. It implies that a pair of cr-particlcs 
can be prepared in two orthogonal states that differ by 
the total topological charge. Computing a product a x 
■ ■ ■ X a for 2n u-particles one can easily get 

^x2n^2»-il + 2"-iV- (2) 
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FIG. 1: Braid group generators. 

Thus if one creates 2n a-particles from the vacuum, there 

is a 2"~-'^-dimensional subspace of states that can be dis- 
tinguished by fusing some pairs of particles together and 
observing a type of the resulting particles. It is used as 
a computational space in the TQC scheme. 

B. The braid group representation 

Recall that exchange statistics of particles living in the 
2+1 dimensional space-time is described by unitary rep- 
resentations of the braid group rather than the symmetric 
group (because the clockwise and the counterclockwise 
exchanges are not equivalent). The braid group Bn with 
n strings can be formally described by generators hj, bj^, 
J = 1, . . . , n — 1 (see Fig. 1) that obey the Yang-Baxter 
relations 

bj bk = bk bj for \j - fc| > 1, 
bj bj+i bj = bj+i bj bj+i, for j = 1, . . . , n - 2. (3) 

The strings can be thought of as world lines of particles, 
whose initial and final positions are chosen on the a;-axis. 



Exchange statistics of fr-particles is described by the 
spinor representation of the braid group [13] 

^ : B2n^ C/(2"). 

It is constructed using the spinor representation of the 
orthogonal group SO{2n). Let us introduce the Pauli 
operators crj, crj, (t| on n qubits, and auxiliary Majorana 
operators Ci , C2 , . . . , C2n defined as 

C2j-l = (Ti • • • (Xi cr|_i (Xi crj (g) Ij+i (g) • • • (X /„, 
C2j = O-J (g) • • • (g) (t|_i (g) crj (g) Ij+i (g • • • (g /„, 

(4) 

where / stands for the one-qubit identity operator and j 
runs from 1 to n. The Majorana operators obey commu- 
tation rules 

Cp Cg + CqCp = 26pg I, cl = Cp, for any p, q. 

Then the spinor representation of the braid group gener- 
ators bi,..., b2n-i is defined as 

>fi{bp) =exp(^-^CpCj,+i) = -^(/-CpCp+i). 



(Wc have omitted the overall phase of ip{bp), since it is 
irrelevant for quantum computation purposes.) 

The Yang-Baxter relations can be easily verified using 
the following identity: 

{Cq, if g ^ {p,p + l}, 
Cp+1, if q = p, 
-Cp, if q=p+l. 

It says that an exchange of adjacent cr-particles is cqiiiv- 
alent to an exchange of the corresponding Majorana op- 
erators (up to a sign). 



C. Topological charge measurements 

The multi-dimensionality of the braid group represen- 
tation accounts for the fact that there is more than one 
way to fuse 2n cr-particles into the vacuum (or into ip- 
particle). A process in which two adjacent cr-particles p 
and p + 1 are fused together and then a type of the re- 
sulting particle (1 or ip) is observed can be described as 
a projective measurement of an observable 

Fp — i Cp Cp-y. 1 . 

The eigenvalues -1-1 and —1 correspond to the resulting 
particle's type 1 and if) respectively. 

A type of a particle that one would obtain by fusing 
together all 2n a-particles is measured by a parity oper- 
ator 

Q = a\®---®al = (-i)" C1C2 • • • C2„. (5) 

Note that Q commute with the action of any braid group 
element, as well as with observables Fp. This is a mani- 
festation of the superselection rules — any local operator 
preserves the total topological charge. Any state \^) of 
2n cr-particles that can be created from the vacuum obeys 
Q\^) = +\^)- Analogously, l^") can be prepared starting 
from a single ^-particle iff Q = — j^"). 

Remark: Strictly speaking, fusion is a process reducing 
the Hilbert space of states, since it replaces two particles 
by one. To simpliiy the notations we describe fusion as a 
projective measurement. This is justified, since a fusion 
can always be followed by an auxiliary fission process in 
which the resulting 1 or ^ particle is split into a pair of 
cr-particles. 



III. TOPOLOGICAL QUANTUM 
COMPUTATION WITH ISING ANYONS 

The goal of this section is to introduce a computa- 
tional model that captures all features of TQC with Ising 
anyons. We will show that any computation within this 
model can be efficiently simulated classically. Finally, we 
describe a natural encoding of a qubit by cr-particles. 
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A. Formal computational model 

To define a formal model we just need to extract its 
constituents from Section II — the computational Hilbert 
space with a fiducial initial state, a set of unitary gates, 
and a set of admissible measurements. 

The computational Hilbert space of n qubits 

jr„ = (c2)®" 

will be represented by 2n cr-particles. The initial state 
|0) = |0) (8) • • • (8) |0) is prepared by preparing pairs of 
a-particles (1,2),..., (2n — 1, 2n) from the vacuum. 

A set of elementary unitary gates includes nearest- 
neighbors exchange operations 

Bp = ip{bp) = exp (--^ Cp Cp+i^ . 

For any p < q define a non-local exchange operation 

=exp(-jcpc,). (6) 

Its conjugated action on the Majorana operators is 

{Cr, if r ^ {p,q}, 
Cq, if r = p, (7) 
-Cp, if r = q. 

One can easily verify that a non-local exchange is a com- 
position of 0(n) nearest-neighbor exchanges, namely for 

any p < q — 2 one has 

Bp,q = -Bg-i • • • -Bp+1 Bp • • • -Bg_i. 

The operations Bp ^ constitute a set of elementary uni- 
tary gates in our model. We shall refer to them as braid 
gates. 

Finally, a set of measurements inchidcs nearest- 
neighbors two-particle fusion processes, i.e., non- 
destructive projective measurements of observables Fp = 
—iCp Cp+i. For any p < q define an observable 

Fp^q — iCpCq. 

Taking into account that Fp^q = Bp+i^q Fp Sp_,_j^, we can 
also measure eigenvalue of any observable Fp,q. Summa- 
rizing, the formal computational model is as follows: 

• The Hilbert space: = (C^)®", 

• The initial state: |0) = |0) (g) • • • (g) |0), 

• Braid gates: Bp ^ = exp (^—jCp c^) , 

• Measurable observables: Fp^q = —icpCq. 

We shall refer to this list as a TQC model. It will be 
assumed throughout this paper that TQC operations are 
implemented perfectly (a storage of quantum states in 
J^n is also assumed to be perfect). 



B. Classical simulation of TQC with Ising anyons 

The fact that any computation in the TQC model 
can be simulated classically follows easily from the 
Gottesman-Knill theorem, see [19]. Indeed, taking into 
account the relation Eq. (4) between the Pauli matrices 
and the Majorana operators, and the conjugated action 
of the braid gates Eq. (7), one can easily prove that any 
braid gate maps Pauli operators to Pauli operators under 
a conjugation. Thus all braid gates belong to the Clifford 
group. Since the set of measurable observables includes 
only Pauli operators, we can directly apply the stabilizer 
formahsm [20, 21] to simulate the TQC. 

Another way to deduce the same result is to relate 
the TQC model and the Fermionic Linear Optics (FLO), 
see [22-24]. A theorem proved in these papers asserts 
that any computation within the FLO model can be ef- 
ficiently simulated classically. In terms of FLO opera- 
tions, the initial state |0) is the Fock vacuum, the braid 
gates Eq. (6) are just special case of Bogolyubov canon- 
ical transformations, while the observables Fp,q measure 
single-mode occupation numbers. Then the classical sim- 
ulatibility of the TQC model follows directly from [24] . 

Loosely speaking, TQC with Ising anyons is an inter- 
section of two computational models known to be classi- 
cally simulatable: the Clifford group/stabilizer formal- 
ism model, and the FLO. This is the reason why we 
need two types of "computationally universal" ancillary 
states. The ancilla [04) takes us beyond the Clifford 
group model, while the ancilla \as) introduces a non- 
linearity necessary to go beyond the FLO model. 

In the remainder of this subsection we explicitly de- 
scribe a set of unitary operators and a set of quantum 
states that can be achieved by TQC operations. 

Let G C J7(2") be a group generated by braid gates 
Bp^g for 2n a-particles. To describe G note that a sub- 
group H c G generated by double exchanges Bp = 
—CpCp+i, p = l,...,2n — 1, coincides with the set of 
all even products of Majorana operators (we do not care 
about the overall phase of operators). Thus if one param- 
eterizes a product of Majorana operators c^^ ■ ■ ■ £3^^" by a 
binary 2n-bit string {xi, . . . , X2n), we get H = (Z2)^"~"^. 
Moreover, the subgroup H is normal: BpHB^ = H for 
any p. One can easily check that the factor group G/H 
coincides with the permutation group S2n of 2n objects. 
Thus G can be represented as a semidirect product: 

G = (Z2)2"-^ K S2n. 

To characterize the set of quantum states that can be 
prepared by TQC operations, note that the initial state 
|0) is a stabilizer state with a stabilizer group 

S = (ciC2, C3C4, . . . , C2„_lC2„). 

Applying any sequence of braid gates Bp^q to this state 
is equivalent to updating the stabilizer group according 
to Eq. (7). A new stabilizer group is 

S' = (Cp(l)Cp(2), Cp(3)Cp(4), . . . , Cp(2„_i)Cp(2„)), (8) 
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where p is a permutation of the numbers {1,2,..., 2n}. 

Let IV') be any state with a stabilizer group S' as above. 
A measurement of an observable Fp^q has non-trivial ef- 
fect on lijj) only if Cp Cg is not a stabilizer of \ip). In this 
case p and q must belong to different pairs, i.e., CrCp and 
CqCa are stabilizers of for some integers r ^ s. More- 
over, these arc the only generators of S' that anticom- 
mute with Fp^g. Therefore, measuring eigenvalue of Fp^g 
is equivalent to updating the stabilizer group according 
to 

(. . . , Cj-Cp^ ^q^si • ' •) ^ (• • • ; ^T-Cg, CpCq^ . . .). 

We conclude that any state one can get from the ini- 
tial state |0) by TQC operations can be described by a 
stabilizer group Eq. (8) for some permutation p G S2n- 

Remark: In the arguments above we have ignored 
eigenvalues associated with stabilizer operators which 
may be either +i or —i. Naturally, after each trans- 
formation one has to update the eigenvalues as well. For 
simplicity we skip these details. 



C. Representation of a qubit 

So far we represented a single qubit by a pair of a- 

particlcs. Although this is the most efficient representa- 
tion in terms of resources, it has some serious drawbacks. 
Since a pair of cr-particles prepared in the basis states |0) 
and |1) has the total topological charge 1 and ip respec- 
tively, a qubit cannot be prepared in a superposition of 
the basis states, e.g., |0) ± because they violate the 
superselection rules. 

For this reason we shall represent a logical qubit by a 
group of four cr-particlcs. The basis states |0) and |i) of a 
logical qubit will be identified with physical states |0, 0) G 
and S ^2- Both these states have the trivial 

total charge. A computational subspace spanned by |0, 0) 
and |1, 1) can be specified by an eigenvalue equation 

-C1C2C3C4 IV') = IV')- (9) 

Logical Pauli operators a^, (t^, and acting on the 
computational subspace can be chosen as 

(T^ = —iciC2, 
= -«C2C3, 

= — 2CiC3. (10) 

Clearly, logical Pauli operators can be implemented by 
braid gates, for example, corresponds to winding the 
particle 1 aroimd the particle 2. Besides, TQC operations 
allow one to measure an eigenvalue of the logical one- 
qubit Pauli operators. 

Note that any four-particle braid gate commutes with 
the parity operator C1C2C3C4, i.e., it implements some log- 
ical one-qubit gate. To find a subgroup of J7(2) generated 
by these gates, it suffices to consider the braid gates Bi,2, 



-82,3, -83^4. In terms of logical Pauli operators one has 

Bi,2 = ^3,4 = exp (-1^ a') = e'*^ ( J i ) ' 

B,,3=expH^a-) = i=(_l. - 

(By abuse of notations, we identify a braid gate and the 
corresponding logical operator.) These gates generate 
the one-qubit Clifford group Cl{\) C C/(2). 

The four-particle qubit representation also has some 
drawbacks which come out if one considers two logical 
qubits. Let us show that any two-qubit logical state 

1^) = a |0, 0) 6 10, 1) c |1, 0) |1, 1) e T^ 

that can be prepared by TQC operations has a product 
form: 

1^) = 1^1) ® IV'2). 

Here and \i\)2) are some logical one-qubit states. In- 
deed, we already know that ji/i) obeys stabilizer equations 

Cp(l)Cp(2) = ±i IV'), . . . ,Cp(7)Cp(8) IV') = ±i ll/i), (11) 

for some permutation p e Sg, see Eq. (8). On the other 
hand, the assumption that jV") is a logical two-qubit state 
implies that 

-C1C2C3C4 IV') = -C5C6C7C8 IV") = IV')- (12) 

Obviously, Eq. (11) and Eq. (12) are consistent with each 
other iff for any 1 < J < 4 one has 

p(2j-l),p(2j)e {1,2,3,4}, 

or 

p(2j-l),p(2j)e {5,6,7,8}. 

In other words, each stabilizer Cp^2j-i)Cp{2j) of the state 
[V^) is composed either from generators Ci, 62, 63, 64, or 
from the generators C5, ce, 67, cg. It means that IV") has a 
product structure |V') = IV'i) IV'2)- 

Summarizing, the four-particle qubit representation al- 
lows one to prepare qubits in a superposition, but no 
entangled states can be prepared topologically. 
No-entanglement rule: The only logical states that 
can be prepared by TQC operations from the initial state 
|0) are products of one-qubit states. 

IV. PURIFICATION OF THE 
EIGHT-PARTICLE ANCILLAS 

One way to get around the no- entanglement rule is to 
use some very noisy non-topological operations to pre- 
pare a state p that approximates some logical entangled 
'target' state. Then one can try to improve accuracy of 
the approximation by running a purification protocol in- 
volving only TQC operations. This is the strategy that 
we shall follow in this section. 
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A. Outline 

A target state which we would hke to purify is the 
maximally entangled two-qubit logical state 

lag) = -^(|0,0) + |1,1)) = -^(|0,0,0,0) + |1, 1,1,1)). 

(13) 

It consists of eight <7-particles [32] The quasiparticles 
1,2,3,4 and 5,6,7,8 represent the first and the second 
logical qubit respectively. 

Let us denote -D(H) the set of all (mixed) quantum 
states on the Hilbert space 71. Let p G 0(^4) be 
eight-particle mixed state that we can prepare by non- 
topological operations. A precision up to which p ap- 
proximates \as) can be characterized by a parameter 

e = 1 - (aslplas). 

It will be referred to as an error rate. 

The only assumption we made about p is that it has a 
support only on the even subspace of i.e., 

Qp^pQ = p, (14) 

where Q is the total parity operator, 

Q = C1C2C3C4C5C6C7C8 = o-j (g) cr| (g) C7| (g) (t|. 

This assumption is justified if p is prepared starting from 
the vacuum by a local operator. As was mentioned in 
Section II, the operator Q measures the total topologi- 
cal charge (1 or tp), so Eq. (14) is a consequence of the 
superselection rules. 

An orthogonal projector onto \as) looks as 

\a8){a8\ = ^{I + S,){I + S2){I + Ss){I + Q), 
where 

51 = -C1C2C5C6, 

52 = -C2C3C6C7, 

53 = -C1C2C3C4. (15) 

Given a binary string s = {31,82,33), Sj G {0,1}, con- 
sider a normalized vector e ^4 such that 

5,- = (-1)^^- I*.), i = 1,2,3, Q|*s) = |*s)- 

Notice that jag) = \Sooo) = I'S'o). Obviously, {|^'s)} con- 
stitute an orthonormal basis of the even subspace of ^4. 
Therefore, p can be written as 

p = ^Pst|*s)(*t|, Poo = l-e. (16) 

s,t 

By analogy with quantum error correcting codes, the op- 
erators Sj and the string s will be referred to as stabilizers 
and a syndrome respectively. 



The goal of a purification is to prepare one copy of 
jag) with an arbitrarily small error rate e' starting from n 
noisy copies of \as) with an error rate e. Performance of a 
purification protocol can be characterized by a threshold 
value of e below which the purification is possible, and 
efficiency, i.e., an asymptotic behavior of n = n(e, e') 
for e' ^ 0. Wc shall describe a protocol for which the 
threshold error rate is 

Ss « 0.384. (17) 

and 

n{e,e')^C{-\oge'f (18) 

for any fixed e < Ss and e' — » 0. Here C is a function of 
e only. The protocol succeeds with a probability at least 
1/2 and there is a flag that tells us when it fails. 
The protocol involves the following steps: 

• Dephasing: make p diagonal in the basis {|^s)}; 

• Syndrome whirling: make the probability distribu- 
tion of the non-zero syndromes s ^ uniform; 

• Puriflcation: convert two noisy copies of jag) into 
one clean copy by postselective measurements on 
four pairs of cr-particles. 

In order to achieve an arbitrarily small error rate, these 
steps have to be repeated sufficiently many times in a 
recursive fashion. Below we describe the protocol on a 
more technical level. 



B. Dephasing 

Let 5 be a group generated by 81,82,83. It consists 
of eight elements. Consider a quantum operation 

ues 

It symmetrizes a state over the group 8, thus implement- 
ing a dephasing in the basis {|\['s)}. The stabifizer op- 
erators Sj themselves can be implemented using braid 
gates, for instance, 5i = — 5^2^5 6- Accordingly, 
can be implemented using braid gates, if C/ G is drawn 
randomly according to the uniform distribution. Obvi- 
ously, for any state p one has 

*s(p) = ^p(s)|*s)(*s|, P(S) = (*s|p|*s). (19) 
s 

We shall assume that each ancilla is acted on by $5 be- 
fore it is fed into the purification protocol. It allows one 
to identify quantum states with probability distributions 
of syndromes. 



7 



C. Syndrome whirling 

A probability distribution of syndromes p{s) can be 
brought by braid gates into the standard bimodal form 



p{s) 



1-e if s= (0,0,0), 
e/7 if (0,0,0). 



(20) 



To achieve this, wc will firstly show how to implement a 
cyclic shift on the set of seven non-zero syndromes s ^ 0. 
Then we shall implement a random cyclic shift. 
Consider a braid-gate 



Ui2 = -82,3 bI 7. 



(21) 



Its conjugated action is as follows (only non-trivial part 
of the action is shown): 
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UI2 = < 



( C2 

C3 

C7 



CS, 
-C2, 
-C7, 
C6. 



Accordingly, a conjugated action of U12 on the stabilizers 



Si - 


-> Si S2, 


S2 - 


- ^2, 


Ss - 


- S3. 



Sj is 



C/12 • ui, -- 

Therefore, U12 implements a XOR-like transformation 

U12 |*si, 82,53) = l*sieS2,S2,S3)5 

where © stands for the addition by modulo 2. Analo- 
gously, one can check that braid gates 

f^23 = -B| 2 -83,4 and U31 = Bi,5 B2,e 
implement XOR-like transformations 



If 



Sl,S2tBS3,S3 



(22) 



?723|4'.i,S2,.3> 

Consider now a braid gate 

U = U3iU23Ui2. 



One can easily check that U implements a cyclic shift of 
non-zero syndromes: 



f/|M/s) = |4',(s)), V 



01234567 
03745621 



Here 77 is a permutation of the numbers {0, 1, . . . , 7}, and 
the syndromes are represented by integers according to 
s = 4si -|- 2s2 + 53- Consider a symmetrization over 
the cyclic group generated by U (it is the cyclic group 
Z7, since 1/"^ = I), i.e.. 



1 ^ 



An application of to a state p = X]sP(^) l^s)(f s| 
transforms the probability distribution p{s) into the stan- 
dard form Eq. (20), where 

e = l-p(O) = 1- (*o|p|*o). 

By construction, $[/ is a probabilistic mixture of braid 
gates. 



D. Elementary purification round 

Suppose we are given a supply of states 

p = ^p(s)|*s)(*s|. 
s 

Here p{s) is some fixed probability distribution of syn- 
dromes (which may or may not have the standard bi- 
modal form). 

Consider a state p p where the first and the sec- 
ond copy is composed from generators ci , . . . , cg and 
C9,...,Ci6 respectively. Let us reshuffle the generators 
by a braid gate B shown on Fig. 2 (inside the dashed 
rectangle) and then measure eigenvalues of four opera- 
tors 



Ti = -ZCgCio, T2 = -iCiiCi2, 
T3 = -iCi3Ci4, T4 = -iCi5Ci6. 



(23) 



Let ti,t2,t3,t4, e {0,1} be the measurement outcomes, 
such that Tj has an eigenvalue (— 1)*^ . Since the in- 
put state p (8) p is a probabilistic mixture of pure states 
l^'r) <?> I^s)) it suffices to analyze the effect of the braid- 
ing-|- measurement operation on those input states. For 
any string of outcomes t — (ti, ^2, ^3, ^4) consider the final 
(unnormalized) state 



where 



p=0 



is the projector corresponding to the outcomes t. Taking 
into account an identity 

(C5C6C7C8)(C9CioCnCi2) = B^ (T1T2T3T4) B (24) 

and the fact that |^'r), |*s) are even states, we conclude 
that 

\Pt\rs) = unless rs ® S3 = ti ®t2 ®t3® U. (25) 
Thus a bit 

t = ti®t2®t3®ti 

can be regarded as a check sum for the syndrome bits 
and S3. If t = then either both syndrome bits are cor- 
rect, r3 = S3 = 0, or both of them are wrong, rs = S3 = 1. 
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tj t2 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 



FIG. 2: The elementary purification round. 

If the input state p has a sufficiently small error rate (the 
probability distribution p(s) is concentrated at s = 0), 
the former possibility is more likely than the latter one. 
As we shall see now, one can enhance a probability of 
states with a correct eigenvalue of ^3 by discarding the 
final state whenever the outcome t = 1 is observed. 

Indeed, suppose we have measured t = and |i^t|rs) 7^ 
0. After some algebra one gets 

\Ft\r.) = \^v)®\tl,t2,t3M) (26) 

(up to a normalization), where 

Ui = ri ® si ® ti ® ^2 ® 1, 

U2 = r2®.S2®t2©i3©l, 

M3 = ?'3 = S3. (27) 

Note that the syndrome bit U3 depends only upon r and 
s, while Ml and U2 depend also upon t. Let us apply 
additional braid gates 

C2C3 : l^'ui.Ma.Ma) ~^ l*Ml®l,M2,M3)> 
C1C2 : |*«i,tl2,«3) ~^ I*tll,ti2el,tl3)- 

conditioned on bits ® ^2 © 1 and ^2 ® ^3 ® 1 respectively. 
These braid-gates flip the bits ui and U2, so we get 

Ml = ri ® si, U2 = r-2 ® S2, Ms = ''a = ss- (28) 

(To avoid clutter, the additional braid gates are not 
shown on Fig. 2.) Summarizing, the output state of the 
elementary purification round is |^'u) with u determined 
byEq. (28). 

For the mixed input state p ® p the syndromes r and 
s are drawn from a product distribution p(r)p(s), so the 
output (normalized) state is 

Pout=^Pout{n)\^^){^J„ 

U 



where 

Pout{^x) = Z-'Y.'^lM^)p{s), (29) 

r,s 

and 

Normalizing pout one gets 

r,s,u r,s 

Note that Z is equal to the probability to observe t = Q, 
i.e., a success probability of the elementary purification 
round. 



E. The protocol 

Let e^-') and e^^^^ be a probability to observe Sj = 1 for 
the distribution p(s) andpou<(s) respectively {j = 1, 2, 3). 
They can be regarded as error rates for the individual 
syndrome bits. If p{s) has the standard bimodal form 
Eq. (20), then e^-^') = 4e/7. On the other, for e < 1 one 
can easily find from Eq. (29) that 

eit\ « 16eV49, e^^^ = e^'j, « 4e/7. 

It tells us that the error rate e^^^ is suppressed quadrati- 

cally, e^'^^j w (e'"'^)^, while the error rates e'^-'^', e^^' remain 

practically unchanged, e^g^^ ~ e'^^\ e^^j « e^^'. For this 
reason we shall iterate the elementary purification round 
shown on Fig. 2 three times to purify all three syndrome 
bits si, S2, S3. The iterations are interlaced with an addi- 
tional braid gate C which shifts the syndrome bits cycli- 
cally, i.e., 

C5i (7^ = 52, (752(7^ = 53, C S^C^ = S^. 

These equations can be satisfied if C transforms the gen- 
erators ci , . . . , cg according to 



ci - 


C6, 


C5 - 


-C7 


C2 - 


C2, 


C6 - 


C3, 


C3 - 


Ci, 


C7 - 


-64 


C4 - 


h, 


Cg - 


Cg. 



An explicit implementation of C is shown on Fig. 3. 



A single round of ag-purification protocol is shown on 
Fig. 4. Its input consists of eight copies of a noisy |ag) 
state in the standard bimodal form with an error rate e: 

p=(l-e)|*o)(*o| + y 

s#0 



1 2 3 4 5 6 7 8 




1 2 3 4 5 6 7 8 



FIG. 3: A braid gate C implementing a cyclic shift of stabi- 
lizers Si ^ §2 ^ ^ Si. 

The protocol outputs a single copy of a noisy jag) state in 
the standard bimodal form with an error rate eout- The 
triangles labeled by 'E' denote the elementary purifica- 
tion rounds shown on Fig. 2. The boxes labeled by 'C 
denote the braid gate shown on Fig. 3. The circle labeled 
by 'W stands for the syndrome whirling transformation. 
Each line on the figure represc;nts c;iglit a-particles. 

The corresponding recursive flow equation eout{^) can 
be found by iterating Eq. (29) three times with an addi- 
tional cyclic shifts inserted after each iteration. Equiva- 
lently, eot,f(e) is implicitly defined by equations 



^out 


= i-z-Vi(o), 


Pi(u) 


= EQ"sP2(r)p2(s), 

r,s 


P2(u) 


= 5]A;^_,p3(r)p3(s), 

r,s 


P3(u) 


r,s 


z 





u 



The initial distribution p(s) has the standard bimodal 
form Eq. (20) with the error rate e. The coefficients 0"g 
and A" 5, are obtained from r"^ by a cyclic shift of in- 
dexes, 

^r,S ~ ^T-L,Ul ^Sl,Ul ^r2ffiS2,«2 ^r-3®S3,U3J 
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FIG. 4: A single round of ag-purification protocol. Time 
flows upwards. Each line represents one copy of the noisy 
lag) state (eight cr-particlcs) . Each triangle E corresponds 
to the elementary purification round shown on Fig. 2. Each 
rectangle C represents a braid gate that shifts the generators 
Si, S2, and S3 cyclically, see Fig. 3. Finally, a circle W stands 
for the syndrome whirling transformation. 



w 0.4 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
e 

FIG. 5: Input /output error rates for a single round of as- 
purification protocol. 

recursively to achieve arbitrarily small error rates. For 
e <C 1 one can easily get 



The final cyclic shift can be discarded because it is fol- 
lowed by the syndrome whirling. We have found a so- 
lution of Eq. (30) using MAPLE. A plot of a function 
eoMt(e) is shown on Fig. 5. 



The threshold error rate 5^ satisfying €out{5s) = 5^ 
turns out to be ~ 0.384. If the initial error rate is 
below the threshold, e < 5%, one can invoke the protocol 



A probability for all elementary purification rounds on 
Fig. 4 to succeed is given by normalizing coefficient Z in 
Eq. (30). For small e one has 

Z = l- d,e + 0{e^). (31) 

The function Z{e) is monotone decreasing on the interval 

< e < ^8 and Z{5s) k 0.04. 

The initial supply of states jag) with an error rate 
eo = e will be called level-0 ancillas. Accordingly, level-fc 
ancillas are obtained from the level-0 ancillas by iterating 
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the protocol shown on Fig. 4 k times. Let and be an 
error rate and the total mimber of level-A; ancillas. The 
numbers e^+i, Uk+i and et, Uk are related by recursive 
flow equations 

Zi^k) 48 2 

rife+i « — g— Uk, ek+i « — e^. (32) 

Here fluctuations of the quantity arc neglected. Monte 
Carlo simulation of the ag -purification protocol shows 
that taking into account fluctuations does not change the 
answer significantly, see Fig. 6. 

If one needs to prepare one copy of \as) with an error 
rate e', the required number of levels k can be found from 
an equation 

~ log (Ceo)' 49- 
The corresponding number of level-0 ancillas is 

fc-i 

no^8'l[z{e,)-\ (33) 

Assuming that k ^ 1 and denoting p(eo) = YljLo ^i^j) 
(one can easily check that this product is convergent), we 
get 



p(eo)log' (Ceo)' 

To find the probability for the protocol to convert no 
copies of the level-0 ancillas into one (or larger number) 
level-fc ancilla, we used numerical simulations, see Fig. 6. 
The success probability Pg = Prob(nfc > 0) was calcu- 
lated as a function of no using Monte Carlo method. For 
each particular k = 1,...,5 an equation Ps(no) = 1/2 
has been solved to find no as a function of eo. As one 
can see from the figure, the scaling of no is pretty well 
described by Eq. (33). 

The operational cost of the purification, i.e., the total 
number of braid gates and fusions needed to achieve an 
error rate e' has the same scaling as no, i.e., it is propor- 
tional to (- loge')^. 



Suppose we have to prepare a large number L of an- 
cillas lag)- Let us first prepare n = 3L • no(eo, e') level-0 
ancillas, split them into 3L groups and then perform the 
og-purification protocol independently in each group. If 
the purification succeeds in each group with a probability 
1/2, the average number of successful group is 3L/2. Us- 
ing the Chernoff bound one can easily show that the prob- 
ability for the number of successful group to be smaller 
than L is at most exp (— L/12). Thus one can say that 
a preparation of a single ancilla jag) with an accuracy e' 
costs about (— log e')^ elementary operations (TQC op- 
erations and preparations of ps)- 




0.05 0.1 0.15 0.2 0.25 



level-0 error rate, eg 

FIG. 6: The graph shows the number of level-0 ancillas, no, 
that one needs to prepare one level-fe ancilla {k — 1, . . . , 5) 
with a success probability 1/2. The success probability has 
been evaluated using Monte Carlo simulation of the protocol 
with 10® trials. Solid lines show the dependence no(eo) that 
one gets using a naive equation no = 8*^ 11^=0 ('^ 
neglects fluctuations of n^). 

This observation also shows that the purification can 
be described by a trace preserving completely positive 
linear map (in the exponentially rare events when the 
purification fails, one can output an arbitrary state). Ac- 
cordingly, we can generalize all above results to the case 
when the preparation of the level-0 ancillas is a stochas- 
tic process that outputs a state pa with a probability Pa, 
such that P'^ Pot = P- 

V. IMPLEMENTATION OF THE CLIFFORD 
GROUP GATES 

Having prepared a supply of clean ancillas jag) one can 
proceed to the next goal — implementation of entangling 
two-qubit gates. We shall now explain how to implement 
a two-qubit controlled gate 

A(a^) :|a,6)^(-l)«''|a,6) 

acting on the logical qubits. Together with logical one- 
qubit Clifford gates which can be implemented by braid 
gates, see Section 111, it will allow us to execute any Clif- 
ford group computation (on the level of logical qubits). 
We shall need the following technical result. 

Lemma 1 The following operations can simulate one 
another with assistance of TQC operations: 

01. A preparation of\a%), 

02. A non- destructive measurement of an observable 
CpCqCrCs (all four labels are distinct). 

03. A unitary gate exp {i^CpCqCrCg) , 
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Remarks: (i) It is meant that one copy of any opera- 
tion can be exactly simulated by one copy of any other 
operation, (ii) The operator CpCqCrCg has eigenvalues 
±1, so 02 can be described by orthogonal projectors 
(1/2) (7 ± CpCqCrCs). (iii) If 01 is not ideal, so that las) 
has an error rate e, then 02 and 03 can be executed 
with an error probability 0(e). (iv) Explicit simulation 
protocols are given in the proof of the lemma. 

The controlled can be easily reduced to 03. Indeed, 
suppose the first qubit is encoded by Ci, . . . , £4, while the 
second qubit is encoded by C5,...,C8. Then A((T^) = 
exp — (jf )(/ — o-f)), where (7| are the logical Pauli 
operators defined as af = —10304 and erf = —ic^ce, see 
Section III. Therefore 

A(cr^) = e** exp (^—i'^c^CiC^ce^ 

■ exp jC3C4^ exp ^csce) . (34) 

The last two exponents in Eq. (34) are braid gates, so 
the controlled cr^ gate is equivalent to 03 (wc disregard 
the overall phase). One remains to prove the lemma. 

Proof of Lemma 1: 
03 can simulate 01: Using solely braid gates one can 
prepare a state with a stabilizer group 

S = {-iciCr, -ic2C8, -10305, -10406) 

(all eigenvalues are +1). Let this state be acted upon by 
an operator 

U = exp(i^ciC2C3C6). 

The conjugated action of U on the Majorana operators 
is as follows: 



-ZC2C3C6, C5 

iCiCsCQ, 06 

-iciC2C6, Or 

C4, Cs 



C5, 

iciC2C3, 

C7, 

Cs- 



Accordingly, the stabilizer group S is mapped into 

S' = {-C2C3C6Cj,CiC3C6Cs, -C1C2C5C6, -C1C2C3C4). 

It coincides with the stabilizer group {81,82,33) of the 
state |os), see Eq. (15). Therefore the state stabilized by 
8' coincides with jas) up to an overall phase. 
02 can simulate Ol: (This part is not necessary for the 
proof, but we shall use this result later.) Using solely 
braid gates one can prepare a state with a stabilizer group 

8 = {-iCiC5, 10206, -10307, 10408). 

Let us measure an eigenvalue of -0506070.8 on this state. 
To find the final stabilizer group, choose generators of S 

as 



After the measurement the last generator is replaced 
by — C5C6C7C8 (may be with the opposite sign), which is 
equivalent to a generator — C1C2C3C4. The resulting stabi- 
lizer group coincides with the one of |as). The prepara- 
tion of |as) by measuring cscecycs is illustrated by Fig. 8. 
01 can simulate 02: Assume that one copy of jos) is 
available. A sequence of braidings and fusions that allows 
one to measure an eigenvalue of C1C2C3C4 is shown on 
Fig. 7. The particles labeled by 1,2,3,4 on the figure 
are prepared in an arbitrary initial state li/'in)- The state 
I as) is described by generators C5, . . . ,612. Accordingly, 
the circuit shown on Fig. 7 is apphed to a state IV'm'Sias)- 
The particles are reshuffled by a braid operator and then 
observables T\ = —iciC2, T2 — —10304, T3 = —ic^ce, and 
T4 = —10708 are measured. The final state |V'/) is read 
out from the particles 9, 10, 11, 12. 

Let (—1)*^ be the measured eigenvalue of Tj. We shall 
consider in details only the case ti ® i2 © ^3 © ^4 = 0. 
Let |$t) be the final state corresponding to outcomes 
t = {ti,t2, t3, 14). Using the stabilizer description of |as), 
see Eq. (15), one can easily check that 

l^t) = (4^ 4^)(4^ 0% Cl\ C*/2)|$o) 

whenever ii ^2 © ^3 © ^4 = 0. Let us apply additional 
braid gates £269, C4C10, cecu, and csCi2 controlled by 
classical bits ti,t2,t3, and t4 respectively (these gates are 
not shown on Fig. 7 to avoid clutter). They map |$t) into 
|$o), so it suffices to analyze the case tj = 0. Obviously, 
|$o) has a product structure |<I>o) = |0, 0,0,0) (g) |V'/m)- 
One can easily notice that the left- upper part of Fig. 7 
with tj = is almost identical to the preparation pro- 
cedure for (as I, see Fig. 8. The only missing element is 
a projector (l/2)(/— cscgcycg). However, one can safely 
add the missing projector because it stabilizes the state 
|as). Therefore, for the outcomes tj — the protocol 
shown on Fig. 7 coincides with the one shown on Fig. 9. 
Taking into account that |as) is the encoded EPR state, 
|as) = 2-V2(|o,0) + |l,i)), the protocol on FIG. 9 is 
just a projection of \tpin) onto the code subspace fol- 
lowed by tclcportation of the encoded qubit from the 
particles 1,2,3,4 to the particles 9,10,11,12. Accord- 
ingly, \^f) = (1/2)(J - ciC2C3C4)\i>in) , up to an over- 
all normalization constant. Using similar arguments one 
can check that \ipf) = (l/2)(/-|-ciC2C3C4)|'^i„) whenever 

ti®t2®t3®t4 = 1. 



02 can simulate 03: (This result has been already 
proved in [30].) Suppose we want to implement an op- 
erator exp (if C1C2C3C4). Let us prepare an ancillary pair 
of particles 5,6 in the state |0). Accordingly, any input 
state IV') of the system satisfies 



8 = {-C1C2C5C6, -02030607,-03040708,10408). 



(C5 + ic6) = 0. 



(35) 
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tj ^2 ^4 



9 10 11 12 




FIG. 7; Implementation of a non-dostructive measurement of 
C1C2C3C4 that consumes one copy of |o8). 
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I a,) = 



^1 6\ 7| 8\ 

(I-C5C(,CjCs)/2 




«> \«> \o> 



FIG. 8: Preparation of log) via eigenvalue measurement of 

C5C6C7C8. 



use the following identities: 

exp (if C1C2C3C4) \ip) = 

= 2exp(fc3C6) n^'^nf IV') 

= 2i exp (f C1C2) exp (§63(54) exp (f cgce) hL'^hL^' 
= 2iexp (f C1C2) exp (IC3C4) exp (-f cgce) U^^^uf 

= 2exp(-fc3C6) nL'^nL'^ IV') 

(36) 

(wc have used Eq. (35)). In each of the four cases one can 
apply a suitable correction operator Uyz (for example, 
U++ = exp(f csce) if the outcomes were ++, etc.) so 
that 

exp (i^CiC2C3C4) IV) = 2[/^,n(2)nW IV;). 

Each of the four outcome combinations occurs with prob- 
ability 1/4. The final state is always the desired one, i.e., 
exp(if C1C2C3C4) IV;)- 



VI. UNIVERSAL QUANTUM COMPUTATION 

The protocols described in Section V allow one to ex- 
ecute any Clifford group gates on the logical qubits. In 
addition to that, one can measure logical qubits in the 
standard basis and prepare fresh logical qubits in the 
state |0). Let us refer to this set of operations as Clifford 
operations. As wc have shown, Clifford operations can be 
implemented with an arbitrarily small error rate and the 
overhead is poly logarithmic. To simplify the discussion, 
we shall firstly set the error rate of Clifford operations to 
zero and then address the precision and overhead issues 
separately. 

Below we will show how to execute the 7r/8-rotation 



Let us measure an eigenvalue of C1C2C4C5. Depending 
upon the outcome, the initial state lip) gets multiplied 

by a projector n^^-* = (l/2)(/± C1C2C4C5) (with a proper 
normalizing coefficient). Next we measure an eigenvalue 
of —16365. The eigenvalues ±1 correspond to projectors 
= (1/2)(1 =F 46365). We claim that after some correc- 
tion depending on the measurements outcomes, the pro- 
tocol effectively executes the operator exp (if 61626364) 
while leaving the ancillary pair of particles intact. The 
correction step requires only braid gates. Indeed, one can 
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FIG. 9: Teleportation. 



A(e' 



1 




e i 



(37) 



on the logical qubit. It is well known that the 7r/8- 
rotation together with Clifford operations constitute a 
universal set of gates. 

Although the 7r/8-rotation cannot be implemented by 
Clifford operations only, we can follow the same strat- 
egy as in Section IV, namely, try to use very noisy non- 
topological operations to prepare a state p that approxi- 
mates some logical target state \a), improve an accuracy 
of the approximation by running a piirification protocol 
(that now can use any Clifford operations) , and then con- 
vert I a) into the gate A(e*T). 

It is not apriori clear what ancillary state |a) leads to 
the most efficient implementation of A(e*^). We shall 
argue that a good choice of \a) is a state 



l«4) 



1 

71 



(|0)+e^^|I))=A(e^^)|+), 



where \+) = 2-'^/'^{\0) + |1)). The state |a4) is com- 
posed from four cr-particles. A purification protocol for 
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|a4) with a high threshold error rate and polylogarith- 
mic overhead which uses only Clifford operation has been 
put forward in [25] under the name "magic states distil- 
lation" . For the sake of completeness wc briefly describe 
it below. Then we assess an efficiency of the whole sim- 
ulation scheme. 

In the rest of this section a word 'qubit' refers to a 
logical qubit encoded by four u-particles as explained in 
Section III. By abuse of notations we shall abbreviate |0) 
to |0) and |T) to |1). Accordingly, I04) will be regarded 
as a one- qubit state. 

A. Converting [04) into a non-Clifford gate 

We start from explaining how to execute the gate 
A(e'T) using Clifford operations and one copy of [04) , 
see [25]. Let [f/') = a\0) + b\l) be an unknown state (the 
coefficients a and b may actually be quantum states as 
well). Suppose we want to apply the gate A(e*5 ) to 
Let us start from a two-qubit state \'^o) = IV' (X" 04) and 
measure an eigenvalue of observable Ti = (S) (recall 
that any multi-qubit Pauli operator can be converted by 
Clifford gates into onc-qubit operator cr^, which is an ad- 
missible observable in the TQC model). The outcomes 
±1 of the measurement appear with the probability 1/2 
each, yielding the final states 

|*+) = a|0,0)+6e't|l,l), 

= ae't|0,l)+6|l,0). (38) 

Applying the controlled operator A(c7^) with the first 
qubit as a control one, we get 

[*+) = A(a-)|*+) = (a|0) + 6e*T|i))^|0), 
1*2-) = A(a-)|*r) = (ae^^|0)+6|l))®|l). 

Now let us measure the second qubit in the {|0),|1)} 
basis. If the outcome is |1), apply additional Clifford gate 
ii" = |0)(0| -I- i |1)(1| to the first qubit (as was mentioned 
in Section III, ii' is a braid gate). In both cases wc end 
up with the final state a|0) -|- 6e*T|l). Thus the input 
state IV') has been acted upon by A(e*7). 

B. Purification of [04) 

Here we outline the magic states distillation method, 
see the original paper [25] for details. A noisy [04) state 
will be described by a onc-qubit dcmsity matrix p. A 
quality of p is characterized by a parameter 

e = 1 - (a4|p|a4) 

which will be referred to as an error rate. The purifi- 
cation protocol exploits some nice properties of the CSS 
second-order punctured Reed-MuUer quantum code. It 
encodes one qubit into 15 qubits and has the minimal 



distance 3. Let 11 be a projector on the code subspace of 
the Reed-MuUer code. Consider a state 

Although po^jf is a 15-qubit state, it can be regarded as 
a one-qubit state encoded by the Reed-MuUer code. It 
turns out that an error rate eout of the state pout is cu- 
bically suppressed as compared to the error rate of p, 

eo«t = 35 e3 + 0(e^). 

The properties of the Reed-MuUer quantiuii c:ode that are 
responsible for this effect are (i) The minimum Hamming 
weight of (T^-type errors that are not detected by the code 
is 3; (ii) The code has non-Clifford automorphisms: an 
operator A(e*^)®^^ commutes with H and its action on 
the encoded qubit coincides with A(e'^). 

This observation provides a natural mean of purify- 
ing p. Namely, one takes 15 copies of p and measures 
eigenvalues of 14 stabilizer operators for the Reed-MuUer 
code. All stabilizers are the Pauli operators, so these 
measurements require only Clifford gates and admissible 
TQC measurements. The final state is accepted iff one 
observes the trivial syndrome (eigenvalue of all stabilizer 
operators is -1-1). After that one applies a decoding trans- 
formation (a certain Clifford group operator) that maps 
Pout into a one-qubit state. The threshold value of e 
is determined by an equation eo„t(e) = e. Denote the 
threshold by ^4. Its numerical value is 

54 « 0.141. 

If e < (54, the output state pout is more clean than the 
input one, i.e., touti^) < e- 

Let Ps be the probability for this algorithm to succeed, 
i.e., the probability to observe the trivial syndrome. In 
the limit e ^ one has Ps ~ 2"^". Moreover, by in- 
troducing an additional "error correction" step into the 
algorithm one can accept a larger set of measured syn- 
dromes (syndromes for which only all cr^-type stabilizers 
have eigenvalue +1). The error correction step enhances 
the success probability to Ps ~ 1 (in the limit e ^ 0). 

The initial supply of states p with an error rate eo = e 
will be called level-0 ancillas. Accordingly, level-fc ancil- 
las are obtained from the level-0 ancillas by iterating the 
elementary purification procedure k times. Let Cfc and 
rik be an error rate and the total number of level-A; an- 
cillas. The numbers Ck+i, rik+i and e/c, rik are related by 
recursive ffow equations 

nfe+i » 77, efe+i«35e^ (39) 

(we are interested in the asymptotic regime e <C 1). Ac- 
cordingly, if one needs to prepare one copy of [04) with 
an error rate e', one needs to have a supply of 

no- |log(e')|^ 7 = logs 15 « 2.5 (40) 

level-0 ancillas with an error rate below the threshold, 
eo < 64. The operational cost of the purification, i.e.. 
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the total number of Clifford gates and standard mea- 
surements needed to achieve an error rate e' has the same 
scaling as no- 



C. Efficiency analysis 

Suppose our goal is to simulate a quantum circuit with 
N one-qubit and two-qubit gates operating on n qubits. 
We assume that the following gate set is used: 

i=(; A(e^f), A(.^). (41) 

The simulation must be able to reproduce the output of 
the circuit with a constant error probability. Accordingly, 
the non-topological gates A(e*^) and A(cr^) have to be 
simulated with an error probability 5 ~ N~^. As we 
have learned in Section IV, preparation of jag) with an 
accuracy 6 requires about (log (5~^))^ raw ancillas ps and 
about the same number of TQC operations. According 
to Section V one copy of jag) can be traded for a gate 
A(cr^) implemented with about the same precision. Thus 
each A(cr^) gate 'costs' 0((log (A''))^) TQC operations 
and raw ancilla preparations. 

Simulation of the gate A(e*^) deserves more careful 
analysis. Consider one round of a4 -purification at the 
level k. It takes as input 15 copies of level-fc ancillas 
[04) with an error rate ek and outputs one copy of [04) 
with an error rate €k+i (sometimes it outputs nothing be- 
cause we use postselection). An implementation of this 
a4-purification round requires 0(1) gates A((j^). To sim- 
ulate each of these gates the ag-purification protocol has 
to be invoked. Obviously, at this point it does not make 
sense to purify lag) ancillas all way down to the error 
rate 6 ~ N~^. Instead, the error rate 0(e|) is sufficient, 
since it still gives the flow equation ek+i = Ce| for 04- 
purification with some constant C. Comparing Eq. (18) 
and Eq. (40) one can see that for a fixed error rate the 
simulation of A((t^) is more demanding in terms of re- 
sources than the simulation of A(e'T). Therefore, we can 
try to use the above observation to improve the efiiciency 
of the whole simulation scheme. 

Indeed, purification of one copy of |ag) with the final 
error rate 0(e|) requires mk ^ (log(efc))^ elementary 
operations. Prom Eq. (39) one gets Ck ~ exp(— cS*^), 
where c > is a constant. Therefore, nik ~ S^'^. The 
total number of a4-purification rounds on the level k is 
Qk ~ rik+i ~ nolS"*^"^, where no is the number of level-0 
ancillas |a4). Thus the total number of elementary op- 
eration needed to generate all level- (fc + 1) ancillas [04) 
is Mfc = nikgk ^ nol5~*'3^*'. Clearly, Mk grows expo- 
nentially with fc, so almost all resources needed to pu- 
rify 1 04) are spent at the highest level of a4-purification. 
Accordingly, the total number of elementary operations 
needed to purify one copy of [04) with the final error rate 



S ~ iV-i is 

d 

Mtot = ^ Mfc « Md « 33'', d « logis K). 

From Eq. (40) with 6 ~ N-^ one gets no ~ (log(A^))^. 
Therefore 

Mt„t~(log(Ar))3. 

We conclude that any gate in the universal gate set 
Eq. (41) 'costs' 0((log (N))^) elementary operations. 



VII. IMPLEMENTATION OF 
NON-TOPOLOGICAL OPERATIONS 

This part of the paper is rather speculative, since we 
know almost nothing about non-topological properties 
of anyons, such as effects of a finite separation between 
particles, non-adiabaticity of the anyonic transport, in- 
teraction between an anyon and a control device, e.t.c. 

Recall that we need non-topological operations to pre- 
pare (may be very noisy) ancillary states \a4) and |ag) 
composed from four and eight cr-particles respectively. 
We shall argue below that a good strategy is to use di- 
rect short-range interaction between anyons. One can 
expect that the amplitude of this interaction decays as 
exp (—I/Ih), where / is a separation between the par- 
ticles and Ih is the magnetic length (for experiments 
with AlGaAs/GaAs heterostructures the magnetic field 
corresponding to = 5/2 is B « 5T, so that Ih = 
(/ic/eB)i/2«io-6cm). 

Remark: Note that any state in the orbit of |a4) or 
|ag) under the action of braid gates is equally acceptable 
as the states |a4) or |ag) themselves. As the number of 
particles increases, the size of the orbit grows, and thus 
the set of acceptable states becomes larger. For exam- 
ple, the orbit of \a4) consists of 12 states (we disregard 
the overall phase), while the orbit of |ag) consists of 240 
states [33]. It suffices to prepare any of these states (we 
have to know which) with a fidelity above the threshold 
one. 



A. How to prepare |a4) 

Let us start from preparation of |a4) since it is much 
easier. The preparation process is illustrated on Fig. 10. 
One starts from the vacuum state, creates two pairs of 
(T-particles and then brings two particles, one from each 
pair, sufficiently close to each other. After that one waits 
for a time r and finally returns the particles to their 
original positions. 



Taking into account that the short-range interaction is 
a local operator, we infer that the total charge of the four 
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has to cancel unwanted interactions. For example, if cr- 
particles 1,2,3,4 are sufficiently close to each other, the 
interaction Hamiltonian looks as 



FIG. 10: A preparation of I04) based on the short-range two- 
particle interaction. 



particles and the total charge of the particles 1, 2 must 
be preserved. Therefore, the interaction can be described 
by a Hamiltonian 

Hint = -iciC2 (g) X -h / (g) y, 

where X and Y are some operators acting on the envi- 
ronment. 

A purpose of the two braid operations preceding the 
interaction on Fig. 10 is to create a state 

10) = B2,s |0, 0) = 2-V2 (|o, 0) + |1, 1)) e J^4. 

Using the qubit representation of Section III one gets 

!</))= 2-V2(|o) + |i)). 

Free evolution under Hint for the time r maps into 
a state 

(|0) (g e^(^+^)^ \-^e) + |1> O e^(-^+^)^ \^e)) . 

Here |^'_e) is the initial state of the environment (one can 
always assume that it is pure). Tracing out the environ- 
ment we end up with a mixed state 

The case p = |a4)(a4| corresponds to r = e™^^. By 

varying the interaction time t we can try to fulfill the 
threshold condition (a4|/9|a4) > 1 — 64 ^ 0.86. This may 
or may not be possible, depending upon particular form 
of X, Y, and |4'b). For example, if X is proportional to 
the identity operator, X = gl, one gets 



r = e 



2igT 



Tuning r such that gr = ±7r/8 we can prepare the desired 
state \ai) (or a state that can be converted to \ai) by a 
braid gate). 



B. How to prepare los) 

Preparation of jag) based on the direct short-range 
interaction between anyons is more tricky because one 



3,k 



<x 



jk - C1C2C3C4 



where Xjk, X, and Y are some operators acting 
on the environment. Recall that \ag) can be pre- 
pared by TQC operations and a non-linear gate W = 
exp (i7r/4 C1C2C3C4), see the first part of the proof of 
Lemma 1 . Free evolution under the Hamiltonian H might 
be used to implement W, provided that one can "turn 
off" the quadratic interactions CjCk ® Xjk- In principle, 
it can be done using a technique analogous to decou- 
pling and refocusing in the Nuclear Magnetic Resonance. 
Indeed, denote F = C1C2, G = C1C3 and consider a Hamil- 
tonian 

H' = ^{H + FHF'f + GHG^ + {FG) H {FG)''). 

One can easily check that 

H' = -C1C2C3C4 ^X + I^Y. 

Now let Ur and U!^ be unitary operators describing evo- 
lution under the Hamiltonians H and H' respectively for 
a time r. If r is sufficiently small, one gets from the 
Trotter expansion 

[/; « C/x . (Ff/jFt) . (Gt/jGt) . [FGU^G^'F'') . 

Therefore one could try to simulate U!^ by 11^-^4 and "con- 
trol pulses" F and G. Obviously, F and G can be im- 
plemented by braid gates (for example, F corresponds 
to winding particle 1 aroimd the particle 2). However, 
before applying any of these braid gates one has to re- 
turn the particle 1, 2, 3, 4 into original well-separated po- 
sitions. After that one can compose the evolutions to 
simulate any desired interaction time. 

The preparation of jag) based on the refocusing may 
fail to provide the necessary precision 5s because it in- 
volves too many noisy operations. So it is more fair to say 
that an additional non-topological operation is needed. 

According to Lemma 1, the state jag) can also be pre- 
pared by TQC operations and a non-destructive measure- 
ment of an observable C1C2C3C4. In other words, we have 
to measure the total topological charge (1 or tp) of four a- 
particles without destroying their pairwise correlations. 
It is very likely that such a measurement can be imple- 
mented using an interferometric device proposed recently 
by Bonderson, Kitaev, and Shtengel [16] (see also [18]) 
to test topological properties of cr-particles. 

The device is based on the Hall bar geometry, see 
Fig. 11, so that the transport of electric charge is gov- 
erned by edge currents on the top and bottom edges of 
the bar. Electrical gates are used to create two constric- 
tions in the region occupied by the FQH electron gas (the 
unshaded region on Fig. 11), so that tr-particles can tun- 
nel between the top and bottom edges through the elec- 
tron gas at either constriction. The parameters of the 
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device are tuned to allow quantum interference between 
the two tunneling paths. The total tunneling current is 
measured through the longitudinal resistance Rxx- Ide- 
ally, such a measurement projects the initial state onto 
an eigenvector of the tunneling current operator. 



C1C2C3C4. Combining the interferometric experiment with 
the standard TQC measurements one can calibrate the 
device to infer an eigenvalue of C1C2C3C4 from the mea- 
surement outcome. 




12 3 4 



FIG. 11: A two-point contact interferometer. Four u-particles 
are trapped at antidots inside the interferometer loop. Elec- 
tric current propagates along the top and bottom edges of 
the FQH electron gas (the unshaded region). Tunneling of 
cr-particles occurs at the constrictions. 




FIG. 12: A braid b describing interference of the two tunneUng 
paths in the two-point contact interferometer. 



Suppose that four antidots arc created inside the in- 
terferometer loop, and exactly one cr-particle is trapped 
at each antidot. Let us label the trapped cr-particles by 
1,2,3,4 and the tunneling cr-particlc by 0. The differ- 
ence between the two tunneling paths corresponds to a 
braid b in which the tunneling particle winds around the 
trapped particles 1,2,3,4, see Fig. 12. Using the braid 
group representation described in Section II one can eas- 
ily find that the action of b is ip{b) = +C1C2C3C4. Thus 
the longitudinal resistance measurement projects the ini- 
tial state of the particles 1, 2, 3, 4 onto an eigenvector of 
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